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Abstract 

The dispersion relations in the real and imaginary parts of the meson self-energies are studied 
to check the consistency of the " renormalization" in cutoff field theory. It is shown that the 
dispersion relations are preserved by the "renormalization" even if the finite cutoff and regulator 
' are introduced in the calculation by hand. 
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5^ ■ 1 Introduction 



In the past two decades, nuclei and nuclear matter have been studied in the framework of quantum 
hadrodynamics (QHD).|§ § The meson mean-field theory of nuclear matter Q has produced suc- 
cessful results to account for the saturation properties at normal nuclear density. Following these 
successes, many studies and modifications have been performed on relativistic nuclear models. One 
of these modifications is the inclusion of vacuum fluctuation effects, which cause divergences in 
physical quantities when they are naively calculated. Chin|| estimated vacuum fluctuation effects 
in the Hartree approximation by using a renormalization procedure, and found that vacuum fluc- 
tuation effects make the incompressibility of nuclear matter smaller and closer to the empirical 
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i value than in the original Walecka model. 

Although the relation between QHD and the underlying fundamental theory, i.e., QCD, is an 
open question, it is natural that QHD is not valid at very high-energy. In this point of view, 
a cutoff or a form factor should be introduced into the theory of QHD. One may introduce the 
cutoff Q or the form factor] 5| to avoid the instability of the meson propagators in the random phase 
approximation (RPA)j^, (7|. Cohenj8| introduced a four-dimensional cutoff into the relativistic 
Hartree calculation and found that the vacuum energy contribution may be somewhat different 
from that in the ordinary renormalization procedures, if the cutoff is not so large. 

One troublesome problem in use of the cutoff is that physical results depend on a value of the 
cutoff and a form of the regulator which are introduced into the theory by hand. It is difficult to 
determine a suitable value of cutoff and a appropriate shape of the regulator phenomenologically. 
For the high-energy physics near the limitation of the theory of the particle physics, Lepage g 
proposed to use the extended effective action (Lagrangian) to eliminate such a dependence on 
cutoff which is by hand introduced into the theory, rather than to search the phenomenologically 
favorable cutoff and regulator (or form factor). This idea is based on the renormalization in the non- 
perturbative renormalization group (NPRG) method. |L0|-|l^]. In the series of papers @, 
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we have studied nuclear matter properties and the vertex corrections in the framework of the cutoff 
field theory, eliminating the cutoff dependence by using the Lepage's method. 

Although the utility of our method of renormalization has been checked phenomenologically 
PH , one may wonder the finite cutoff can be introduced consistently and does not break the 
consistent condition. In this paper, we study the dispersion relation between real and imaginary 
parts of the meson self-energy in the framework of the finite cutoff field theory to check the 
consistency of our method and show how the consistency is preserved in our formalism. 

This paper is organized as follows. In §2, we review and reformulate our method to remove 
the cutoff dependences of the physical results. In §3, we study the dispersion relation between real 
and imaginary parts of the meson self-energy in the framework of the finite cutoff field theory to 
check the consistency of our method. In §4, the influence of the finite cutoff effects to the vacuum 
polarization is discussed. Section 5 is devoted to a summary. 



2 Vacuum polarization effects in cutoff field theory 

In this section we review the "renormalizations" based on cutoff field theory, discussing the vacuum 
polarization effects. After that, we reformulate the method in the form of differential equations. 

At first, we start from " renormalizable " Lagrangian which defined at cutoff Ao. The cutoff 
Ao denotes a upper limit energy scale below which the theory can be used. We will show that 
"nonrenormalizable" terms will appear in the low-energy effective Lagrangian, even if we start 
with a "renormalizable" Lagrangian. We demonstrate the Lagrangian does not have to be "renor- 
malizable" in the sense of the effective action and should have the " non-renormalizable" terms at 
low-energy. 

The "renormalizable" Lagrangian of the a-uj model is 
L Aa = i>{i^d» - M + g s <f> - g^Vty + \d^d^ - ^m 2 J 2 + ^S^^ - ±F <f 
-U{<j>) - \f^F^ + ImlVnV, 

4 

u(</>) = £)c„G^) B , (i) 

n=0 

where ip, (f>, Va, M, m s , m v , g s and g v are the nucleon field, cr-meson field, w-meson field, nucleon 
mass, cr-meson mass, w-meson mass, c-nucleon coupling and w-nucleon coupling, respectively. The 
vector field strength is F M „ = d^Vv — d u V^, and C n are "bare" cr-meson self-interaction couplings 
which are adjusted to reproduce renormalization conditions. For simplicity, in Eq. ([!]) we omitted 
the other counter-terms which are not necessary for the relativistic Hartree approximation(RHA) 
and the random phase approximation (RPA) for cr-meson self-energy. The Lagrangian ([!]) is valid 
only in the low-energy region below Ao . 

For nucleon propagations, we use the solution in RHA. In this approximation the meson fields 
are replaced by their expectation values, 

(j> > ((/))= 0o 

^ — » (V^)=6^V . (2) 
At zero-density, the nucleon propagator G{p) is just the same as the Feynman propagator for free 
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nudeon, i.e., 

<*■» = ^& P) 

Let us estimate one- loop vacuum polarization in the RPA, using this nucleon propagator. For 
simplicity of the illustration, in this section, we consider vacuum polarization at zero baryon 
density. The extension to the case of the finite baryon density is straightforward as in the ordinary 
renormalization procedure. || We will discuss the RPA calculations at finite density in the next 
section. 

If we use the sharp regulator with the cutoff A , one-loop RPA contribution of the cr-meson 
self-energy (Fig. |l|) is given by 

nV,A ) = ^fdx [ A ° dk' k ' 3 }- k ' 2 + (4) 



2tt 2 7 Jo (k' 2 +A 2 ) 2 
where A 2 = M 2 - q 2 x(l - x), k' = k - q(l - x). 



Fig.l 



In Eq. (^), the cutoff Ao appears. Suppose that we do not know a correct value of Ao which 
means a limit of an energy scale of the theory and we use another cutoff A which is a smaller than 
Ao in the calculation. (Of course, A should be larger than the energy scale of the physics in which 
we are interested.) 

As in the renormalization group equations, we require that physical quantities do not depend 
on A, although the physical quantities depend on Ao. To achieve this, we estimate the contribution 
which is needlessly discarded by using the cutoff A instead of Ao. In the case of cr-meson self-energy 
d), this is given by 

An( g 2 ,Ao,A) = n( 9 2 ,A )-II(g 2 ,A) 
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2ir 2 J a J A (k' 2 + A 2 ) 

<a 2 [\ x [ a dv v i- y+z -r {1 -:l (5) 



4tt 2 Jo Ji [v + z ~ ux(l — x)] 2 

where y = k' 2 /A 2 , z — M 2 /A 2 , u — q 2 / A 2 and a = A 2 , /A 2 . If external momentum q is smaller 
than A, we could expand the Eq. (||) around u — and we get 

00 1 d n 

An(g 2 ,A ,A) = ^-_An| u=oW " 

n=0 



^^A 2 (q 2 /A 2 )". 

oo 

E^«(3 2 )" ( 6 ) 

n=0 
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Each term of expansion (||) corresponds to a quantity which is proportional to (d 2 ) n (/) 2 in La- 
grangian. Therefore, if we add terms, which is equivalent to (^|), to Lagrangian at the beginning 
of our calculation and determine them phcnomcnologically, wc do not need to know about Aq. In 
that way, we get proper physical results by using the cutoff A instead of Ao. 

In other word, this means that the A-dependence of the physical results is removed by the phe- 
nomenological determinations of the new terms, since A is introduced by hand and has no physical 
meaning. To say exactly, we must determine all coefficients of expansion @ phenomenologically, 
to eliminate all A-dependence. It is not possible actually. However, if A is the same order as Ao and 
if q 2 is smaller than A 2 , it is easily seen that D* is of order 1 and the n-th term in the expansion 
) is order of (1/A 2 )™ -1 . Therefore, if we consider a limit of A — > oo, the terms D n (q 2 ) (n > 2) in 
) amounts to zero, and we need to determine only the coefficients D n (n < 1) phenomenologically. 
In this limit, it is clear that this procedure is essentially the same as the usual renormalization 
procedure in which the coefficient Fq is chosen to renormalize the cr-meson mass and the coefficient 
Fi is chosen to renormalize the cr-meson wave function. We must determine Fq + Dq and F% + D\ 
instead of Fo and F\, if we use A instead of Ao- However, the number of the phenomenological 
inputs is unchanged, since we only need to know the sum Fi + Di and do not need to know Fi and 
-Dj, respectively. In this meaning, the effects of quantum fluctuations are "renormalizable" if we 
re-define the counter-terms as F£ = F + D and F* = F± + D\. 

However, since the cutoff Ao may be several GeV in QHD, we can not take a limit A — > oo. In 
this case, the errors of order (I/A) 2 may not be negligible. 

According to Lepage's proposal, we determine not only the coefficient C n + D n (n < 1) but also 
the higher coefficients D n (n > 1) phenomenologically to remove the errors which occurs from the 
finitcness of A. As an example, we determine the coefficients £>2 phenomenologically in addition to 
F * and F*, if we want results with 0(1/A 4 ) errors. This is equivalent to add the terms, which are 
proportional to (d 2 ) 2 (f> 2 , to the Lagrangian. Moreover, if we want the errors which are 0((1/ A 2 ) N ), 
we should use the effective Lagrangian 

N 1 

L A ^L Ao + Y,^D n (d 2 )^ 2 . (7) 

with the cutoff A. Similarly, adding higher terms to the Lagrangian and determining them phe- 
nomenologically, the A-dependence of the physical results are removed order by order. 

It should also be emphasized that the low-energy effective Lagrangian such as (Q) has "non- 
renormalizable" terms, although the original Lagrangian (Q) is "renormalizable". Of course, they 
cause no problem, because of the coefficients of " nonrenormalizable" terms contain a suppres- 
sion factor (I/A 2 )' 1-1 . If we use the "nonrenormalizable" Lagrangian (Q) with A instead of the 
"renormalizable " Lagrangian (|f|) with the cutoff Ao, we get the proper results. 

We also remark that the Ao-dependence of the effective Lagrangian is almostly included in the 
terms which are determined phenomenologically. Therefore, we do not need detailed knowledge of 
cutoff Aq. 

In the discussions above, we assume the Lagrangian is "renormalizable" at Ao. However, the 
discussions above indicate that the Lagrangian may have following " nonrenormalizable terms" even 
at A . 

oo 

4 =i Ao + E/«(W (8) 

where F n is order of (I/A 2 ,)™ -1 . In this case, we need to determine only F n + D n (n < N) and 
put F n + D n = (n > N) if we want results with errors of order (l/A 2 ) N . The number of the 
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phcnomcnological inputs is unchanged as before. Therefore, in general, we consider the extended 
Lagrangian (||) at the initial. 

Next we reformulate the results above by using a differential equation in the finite cutoff theory. 
PH , [l6| In this formulation, it is more clear that the physical results hardly depend on the regulator 
which was introduced by hand. 

Here we consider cr-meson self-energy (|]) again, using the cutoff A instead of the true one 
Ao. The A-dependence of H(q 2 , A 2 ) is determined by the high-energy behaviour of the integrand. 
Therefore, we can easily see that II has the A-dependence of order A 2 by power counting. The 
ambiguity of 0(A 2 ) is too large, since A is larger than the physical energy scale of the problem. 
Therefore, instead of II itself, we consider a second derivative 

n(2>(92 ' A)= (W n(92 ' A) ' (9) 

of it. We can know that Eq. (||) is 0(1/A 2 ) by power counting. If A is sufficiently large, this means 
that n^ 2 ) hardly depends on A. Especially, in the the usual renormalization procedure (A — > oo), 
n^ 2 ) is independent of A. 

To get the self-energy II itself, we regard Eq. (||) as a differential equation and integrate it. 
00 We get 

n 2 (,7 2 ,A) = F>+F[q 2 + dtj^ dvU^(v,A), (10) 

where Co and C\ are integration constants, and c denotes a "renormalization point". If Co and C\ 
are determined phenomenologically, the A-dependence of II 2 is order of 1/A 2 , and if A is sufficiently 
large, we can ignore the dependence. But we may not be able to ignore A-dependence of II2, if A 
is not so large. In this case, we need to repeat the differentiations until the A-dependence becomes 
sufficiently small. If we differentiate II (TV + l)-times with respect to q 2 , we get 

d N+1 

rt N +V(q 2 ,A) = (dg2)Ar+i n(g 2 ,A) (7V + 1 > 2). (11) 

The A-dependence of Eq. ( pT| ) is order of (l/A 2 ) N . We integrate Eq. ( pi] ) as in the case of Eq. 
d). We get 

JV 

IW 9 2 ,A) = Y, F n(l 2 T 

+ f dsi f 1 ds 2 ■ ■ ■ ■ f N ds N+1 rt N+1 Hs N+1 ,A). (12) 

•/ c J c J c 

If Fq ~ F' N are determined phenomenologically, the A-dependence of IIjv+i is order of (1/A 2 ) N . 

Therefore, if A-dependence of physical quantities are not small enough to be ignored, we should 
differentiate the physical quantities by some external parameter until the A-dependence becomes 
small enough. It is clear that this formulation is equivalent to the one proposed by Lepage. Fur- 
thermore, in this reformulation, we easily see that our formulation does not depend on the details 
of the regulator, since we have not assume its explict form. This means that the phenomenological 
determinations of the coefficients of the effective Lagrangian remove the dependence on the form 
of the regulator which is introduced by hand, as well as the A-dependence. 
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3 Dispersion relations in a-u self-energies 

In this section, we ascertain dispersion relations of RPA a-u> self-energies in finite density system. 

We study the RPA meson self-energies in finite density system. The cr-meson self-energy, the 
w-meson self-energy, the a-u) mixed polarization part self-energy are given by 

- -igl J -^Tr[G(k)G(k + q)], 



IW<?) = -ig'tj -^Tr[^G(k) 7v G(k + q)}, 



, , d 4 k 



nf(g) = i9s9vJ-0^Tr[ lfl G(k)G(k + q )], (13) 

where G is nucleon propagator in RHA and is given by 

G(k) = (Y% + M*)[ — 2 — 2 — + ^S(k* Q E*(k))9(k F - |fc|) 

= G F (k)+G D (k), (14) 

where k F is the Fermi momentum. The M* is the effective nucleon mass which is related with the 
vacuum expectation value 4>o of cr-meson field as 

M* = M — . 9s O - (15) 
The effective three dimensional momentum and the effective energy of the nucleon are given by 

= (k°-g v V°,k), 
E*(k) = Vk 2 + M* 2 . (16) 

In Eq. ([Ti]), G F describes the propagation of nucleon and antinucleon. The Gd describes the 
propagation of holes in the Fermi sea and it also include the effects of the Pauli exclusion principle 
for the particle and antiparticle excitations below the Fermi surface. 

Here we substitute Eq. (|l4|) into Eq. (|l3|), we can divide meson self-energy into a vacuum 
fluctuation part and a density dependent part. 

U(q)=IL V (q)+U D (q). (17) 

The density dependent part and vacuum fluctuation part are given by 

nf(g) = -ig 2 s J -0^Tr[G D (k)G D (k + q ) + G D (k)G F (k + q ) 

+G F (k)G D (k + q)}, 
KM = -i9lj ^Tr[ lf ,G D (k) lu G D (k + q ) 

+l^G D {k) lv G F (k + q )+ jf,G F {k) lu G D {k + q)], 



<(<?) = ig.g v I ra^Gi /,)(/( /, +,)]. us) 



G 



and 



UY(q) = -ig 2 s j -0-Tr[G F {k)G F {k + q)], 



*v J ^ 

= {*f--9r)Kh)- (19) 

We easily see that the vacuum fluctuation part of the cr-meson self-energy is equal to the one 
at zero-density (Eq. (^)) if we use effective nucleon mass M* instead of M. This is true for the 
vacuum fluctuation part of the w-meson self-energy. Therefore, as well as in the case of zero-baryon 
density, the vacuum fluctuation part II is divergent when it is calculated naively. So we introduce 
the cutoff A to regularize the divergences and perform the "renormalization" procedure described 
in the previous section, i.e., in vacuum fluctuation part, we examine a dispersion relations using 
the derivative of IT y rather than II^ itself. In this section, we adopt a limit of A — > oo. The effects 
of the fmitencss of the cutoff to the dispersion relation will be discussed in the next section. 

First we consider dispersion relations in the density part H D . On an energy complex plane, it 
is given by 

rD/ . _ 2P f 00 Imn D (uj,q)-uj 



ReU^iq) = — / du, (20) 



where "P" denotes the principal value. The derivation of Eq. (g0|) is presented in Appendix. 

In Fig. ||, we show dispersion relations of density part Eq. (|o|). In this figure, the results 
calculated directly by using Eq. ([l8]) are shown by solid, dotted, dashed, dashed-dotted lines, 
while the corresponding results calculated by using the dispersion relations Eq. ( pp| ) are shown by 
filled circles. We have put the absolute value of the three dimensional momentum q — lQQMeV 
in the calculation. We see that each line agrees with filled circles. The dispersion relations are 
preserved in the density part. 



Fig. 2 



We assume the same relation in the vacuum fluctuation part 11^ . However 11^ has strong cutoff 
dependence and diverges in the limit A — > oo. Samely, the large to contributions do not vanish 
in the w-integral in the dispersion relation and the w-integral diverges. So we can not apply the 
dispersion relations to H v just as it is. Therefore we consider the dispersion relation in the second 
derivative of H v rather than H v itself. 



d 



2 



[Ren V (q)] = =£- / T dcu . (21) 



2P f°° ImU v {uj,q) -lu 



7T 



The numerical calculations for the dispersion relation Eq. (21) in the vacuum fluctuation part 
are shown in Fig. ||. We see that Eq. ( pl| ) is preserved as well as Eq. (p0|). The dispersion relation 
is preserved in the conversing second derivative of the diverging quantity. 
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Fig 



.3 



Here we add a comment which concerns the result of Fig. ||. When we calculated H v using Eq. 
jl9| ) , we have used the Wick rotation. Similar as the dispersion relation, one condition for the use 
of the Wick rotation is that the integrand of the energy integral in U v should not contribute in the 
limit where a complex variable go to oo. However, this condition is not satisfied and H v diverges. 
This means that, to say exactly, we cannot directly apply the Wick rotation to H v . However, we 
know that we can get the justifiable physical results by using the renormalization procedure after 
the Wick rotation. In other words, the ambiguities which have appeared in the use of the Wick 
rotation is canceled by the renormalization procedure. However, the results above indicate that, 
alternatively, we can consider the conversing derivative of the diverging quantity before we use the 
Wick rotation. In fact, this is more clear interpretation for the Wick rotation. Furthermore, as is 
shown in Appendix, the Wick rotation in fco-plane can be consistent with the dispersion relation 
in w-plane. We can interpret the Wick rotation and the dispersion relation on an equal footing, 
by considering the conversing derivatives of the diversing quantities. 



4 The effects of the finite cutoff to the vacuum polarization 

Next, we study the cutoff dependence in a and u> self-energies numerically, and we examine the 
effects of the fmiteness of the cutoff in the dispersion relations Sec. 3. After that, we discusses the 
renormalization in QHD which is defined by the cutoff which is a order of several GeV. 

We examine the finite cutoff dependence in the derivatives of H v numerically. Consider second 
and third derivatives of the vacuum parts of the meson self-energies with respect to q 2 . We call 
them n y ( 2 ) and n y ( 3 ) respectively. The A-dependence of H v ^ are order of (1/A 2 ), while the 
A-dependence of n^( 3 ) are order of (1/A 4 ). Naturally, in the ordinary renormalization procedure, 
the limit, A — > oo, is taken, and the cutoff dependences of H v ^ and n v ( 3 ) vanish. In the case of 
the finite cutoff, we should check an utility of a method of renormalization in Sec. 2 by numerical 
calculations. 

In Figs. § and |, we show U v ^ and n y(3) calculated with several cutoff A = 2, 3, 5 GeV, 
respectively. Naturally, the cutoff A gives the limit of integration area of internal momentum in 
RPA Feynman diagram (Fig. nj). In these results, we see that, although A-dependence of H v ^ is 
large, the A-dependence of n^ 3 ) is sufficiently small and can be neglected. 



Fig. 4 



Fig. 5 



We also show the results with a finite cutoff A w which is the upper limit of the integration in 
the integral Eq. (p0|) . The results show the similar features as in Figs. and||. These results is due 
to that the integrand in n y ( 3 ) decreases more rapidly than one in n^ 2 ) as w-beocomes large and 



8 



the value of H^ 3 ' is hardly affected by the change of the cutoff. In other word, the ^-integration 
in IT^ 3 ) is mainly dominated by the contributions at small uj which is much smaller than A w . 



Fig. 6 



Fig. 7 



Those results indicate that, in QHD with the finite cutoff of order GeV, if we want to the 
physical results which hardly depend on the cutoff, we should remove the errors of order (1/A 2 ). 

For example, in the vacuum fluctuation part of RPA cr-meson self-energy, we expand around 
renormalization point q 2 

nr(g,A) 



n 



M* = M, 



M*=M 



dM* 



M* = M 



an: 



dq 2 

9 2 nr 



dM*dq 2 

i d 3 nY 



3! dM* 3 

i d 4 uY 



. (I 2 - m 2 s ) -f - 
(M 

*=M 
i =m 2 

(Af 



(M* - M) 

is 2 nr 



M* = M 
„2 = „ 



4! dM* 4 



M* = M 



(M* 



2dM* 2 
M){q 2 ->, 

Mf 

Mf 



(M* - M) 2 

[ 

i d 2 nY 



2 (dq 2 ) 2 



i 9 3 nr 



2 dM* 2 c 



M* = M 
q 2 = n 



(q 2 -m 2 ) 2 

Vf 
2 

(M*-M) 2 (q 2 -m 2 s ) 



(22) 



Here, in the limit A — > oo, the first four terms of the right hand side (r. h. s.) in Eq. ( p2| ) 
are diverging. In the ordinary renormalization procedure, we must determine only these divergent 
terms phenomenologically, and we can calculate the remaining terms. However, since the cutoff is 
finite in QHD, the cutoff dependence of the next five terms in r. h. s. of Eq. (22) is not small. The 
figures |6| and indicates that we can not ignore the error of 0(1/A 2 ) in those five terms. Therefore, 
we should determine 9 terms which has the errors greater than 0(1/A 2 ) phenomenologically. 

Finally, in Fig. [| we show the cutoff dependence of fl^y^) from which the first four terms 
in Eq. ( ]2^ ) are subtracted. In Fig. |9|, we also show the cutoff dependence of ^oti/A*) f rom which 
the first nine terms in Eq.(|2^) are substracted. As expected, we see again that we can not ignore 
the errors of 0(1/A 2 ). In QHD, we must determine much more terms phenomenologically than in 
the usual renormalization procedure. 



Fig. 8 



Fig. 9 
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5 Summary 



The results obtained in this paper are summarized as follows. 

(1) We formulated the method of renormalization in the vacuum polarization effects in QHD 
as the finite cutoff field theory. We can remove the errors which arise from the finiteness and 
ambiguity of the cutoff, order by order, by determining the coefficients of the effective Lagrangian 
phcnomenologically. 

(2) We reformulated the method by using a differential equation in the finite cutoff field theory. 
In this reformulation, it is more clear that the physical result hardly depends on the cutoff A and the 
details of its regulator which are introduced ad hoc. In this method, the A-dependence of physical 
quantity becomes small, by differentiating the physical quantity which has large A-depcndcncc with 
respect to the external parameters. After that, we integrate the obtained differential equation by 
the external parameter. The cutoff dependence is removed by the phenomenological determinations 
of the integration constants which appear in the integration. 

(3) We ascertained dispersion relations of RPA a and w self-energies at finite density. In 
vacuum fluctuation part which diverges in the limit A — > oo, we examined a dispersion relations 
in the conversing derivatives of H v rather than H v itself. It is shown that such a derivatives well 
preserves the dispersion relations, although itself is diverging. 

(4) We examined the dependence on the finite cutoff in the meson self-energies numerically, and 
discussed an utility of our renormalization method in QHD. We found that the errors of 0(1/A 2 ) 
which arise from the finiteness of cutoff should not be ignored in QHD where A is order GeV. 
This result indicates that we must determine much more terms in the effective QHD Lagrangian 
phenomenologically than in the usual renormalization procedure. 

It is very interesting to calculate the effective meson masses at finite baryon density by using 
our method of renormalization. It is in progress. 



The authors would like to thank T. Kohmura, T. Suzuki, M. Yahiro, K. Harada, K.-I. Aoki, 
H. Yoneyama and T. Mitsumori for useful discussions and suggestions. This work is supported in 
part by the Japanese Scientific Foundation (Grant No. 09740208). 



Here we derive the dispersion relation in the case of the particle-hole part of the cr-meson 
self-energy. We can derive the relation for the other parts in similar way. 
The cr-meson self-energy in the RPA calculations is given by 



For this purpose, it is convenient to use the particle-hole-antiparticle decomposition [|f^-(l9| of the 
nucleon propagator rather than the Feynman-density decomposition. 
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Appendix 




(23) 



G(k) 




rl-6(k F -\k\) 6(k F -\k\) 1 



Ik* - E*(k) +ie k*-E*(k)~ie k*+E*{k)-ie 
G p (k) + G h (k)+G a (k), 



(24) 
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where G pi Gh and G a represent the particle, the hole and the antiparticle propagations, respec- 
tively. The particle-hole part of the cr-meson self-energy is given by 

nf(?) = -ig 2 s j ^Tr[G p (k)G h (k + q) + G h (k)G p (k + q)] 

= U + (q)+IL_(q), (25) 



where 

n+(g) 



. 2 M ' 



2 f d A k M* 2 + k-(k + q) 



r) 4 E*(k)E*(k + q) 

x{e(k F -\k\)[i-e(k F -\k + q \)]-. 



[k - E*{k) - ie] [k + q a - E*(k + q)+ ie] 
gl f d 3 k 9(k F - \k\)[l - 9(k F - \k + q\)} 



2 J (2tt 



E*(k)E*(k + q) 



x \ 2(E*(k) + E*(k + q)) + {AM* 2 - q 2 ) 



and 

n-(«) 



S*(fc) + g -S*(fc + g)+ie 
d 4 fc M* 2 + k ■ (k + q) 



(26) 



2 /* d 4 * 



r) 4 E*(k)E*(k + q) 

xlB{k F - \k + q\)[l -6(k F - \k\)]- l — — - — l — - 

\ ' Ul 1 ' Ul [k Q + qo-E*{k + q)-ie} [k Q - E* (k) + ie] 

g 2 s f d 3 k 0(k F - \k\)[l -6{k F - \k + q\)\ 



2 J (2tt 



E*(k)E*(k + q) 



x\ 2(E*(k) +E*(k + q)) + (AM* 2 - g 2 



E*(k)-q -E*(k + q)+ie 



Using the relation, 



1 — = P- =F iir5(x) 



x±ie x 



(27) 



(28) 



we get 



ImJA+(q) 



-g^(AM* 2 -q 2 ) 



d 3 k e(k F -\k\)[l-6(k F -\k + q\)} 



(2tt) 3 

xS(E*(k)+q Q - E*{k + q)) 



E*(k)E*(k + q) 



(29) 



and 



ImSL-{ q ) = - g f-{AM* 2 -q 2 )J^- 



d 3 k 0(k F -\k\)[l-6(k F -\k + q\ 



) 3 

x5(E*(k)-q Q -E*(k + q)). 



E*(k)E*(k + q) 



(30) 
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From these equations, we see that ImU-(q , q) = 7mII + (— g , q). 

If |k + q| < |k|, there is no contribution for the /c-intcgral in Eq. (|2S|). On the other hand, if 
|k + q| > |k|, E*{k) ~E*(k + q) < 0. Therefore, we get 

7mII + (g)=0 for q < 0. (31) 

Similarly, we get 

7mII_ (g) = for q Q > 0. (32) 

Since II + has no pole in the upper half-plane of the complex go-plane, we get the dispersion 
relation 

ReU + (q) = Z f (33) 

Similarly, since II _ has no pole in the lower half-plane of the complex go-plane, we get the dispersion 
relation 

7mII_(£j, q) 



JMI_(g) = -- f 
71" J-, 



?0 



-duo. (34) 



Therefore, we get 

RelLf = ReTl + (q) + Rell-{q) 



P f°° ImU+(uj,q) , P f°° 7mkL(u,g) , 
' dw / duj 



P /-°° Jmn+^.g) ^ Pf° ImU + (-uj,q) d ^ 



TT J 0J-q 7T J_ QO W - 90 

P /- 00 7mn+(w,g) , P fM+kq), 
' acj / dw 



ir J uj - q tt J -oj - q a 

2P f°° ImW s h (aj,q)-to 



du), (35) 

tt Jo ^ - % 

The dispersion relations for the other parts of II can be also shown in similar way. For example, 
if we use the Feyman propagator Gp{k) in II, we get the dispersion relation for II y . In that case, 
we easily see that the condition for the dispersion relation in go-plane can be consistent with the 
condition for the Wick rotation in fco-plane, since we decompose 11^ into 11+ and HY_ to apply the 
dispersion relation. 
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Figure 1: Feynman diagram for RPA meson self-energy. The solid line denotes the nucleon and 
the dotted lines denote the meson. 




Figure 2: Dispersion relation for the density part. The solid, dotted, dashed and dashed-dotted 
lines are iMlf , i?eIl£,iMl£ and Reltff , where Il£ = IT^ - lift and Iff = Ilf 2 = Ug as 
q = (q, 0,0). The corresponding results calculated by the dispersion relation Eq. j20| ) are shown 
by filled circles. We put q = 100 MeV in these calculation. 
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Figure 3: Dispersion relation for the second derivatives of the vacuum fluctuation part. The 
solid, dashed lines are the second derivatives of Rel\^ and ReU^ , respectively. The corresponding 
results calculated by Eq. (^) are shown by the filled circles. 
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Figure 4: The A-dependence of the second derivative of U v . The dashed, dotted, dashed-dotted 
and solid line are the results with A = 2 GeV, A = 3 GeV, A = 5 GeV and A = oo, respectively, 
(a) The cr-meson self-energy, (b) The w-meson self-energy. 
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Figure 5: The A-dependence of the third derivatives of H v . The meaning of each line is same as 
in Fig. U (a) The er-meson self-energy, (b) The w-meson self-energy. 
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Figure 6: The A w dependence of the second derivative of II . The meaning of each line is same 
as in Fig. [|. (a) The a-meson self-energy, (b) The w-meson self-energy. 
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Figure 7: The A w dependence of the third derivative of H v . The meaning of each line is same as 
in Fig. H (a) The er-meson self-energy, (b) The w-meson self-energy. 
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Figure 8: The A-dependence of ^^n/\2y The meaning of each line is same as in Fig. |. (a) The 
cr-meson self-energy, (b) The w-meson self-energy. 
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Figure 9: The A-dependence of m*v The meaning of each line is same as in Fig. |. (a) The 
cr-meson self-energy, (b) The w-meson self-energy. 
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